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1. INTRODUCTION 

The problem of initial singularity still remains at the center of modern day cosmology. 
Though Big Bang theory is deep rooted among the scientists dealing with early day cos- 
mology, it is natural to look back if one can model a Universe free from initial singularity. 
An attempt was made by us for several years. In doing so, we studied the nonlinear spinor 
field in a Bianchi-type I universe (BI) in a series of papers [J, Q, 0, S IE S 0, H 1^ • 

Why nonlinear spinor field? It is well-known that the nonlinear generalization of clas- 
sical field theory remains one of the possible ways to over come the difficulties of the 
theory which considers elementary particles as mathematical points. In this approach 
elementary particles are modeled by regular (soliton-like) solutions of corresponding non- 
linear equations. The gravitational field equation is nonlinear by nature and the field 
itself is universal and unscreenable. These properties lead to definite physical interest 
for the proper gravitational field to be considered. We prefer spinor field to scalar or 
electromagnetic ones, as the spinor field is the most sensitive to the proper gravitational 
one. 

Why anisotropic universe? Though spatially homogeneous and isotropic Friedmann- 
Robertson- Walker (FRW) models are widely considered as good approximation of the 
present and early stages of the universe, the large scale matter distribution in the ob- 
servable universe, largely manifested in the form of discrete structures, does not exhibit 
homogeneity of a higher order. In contrast, the cosmic background radiation, which is 
significant in the microwave region, is extremely homogeneous, however, recent space in- 
vestigations detect anisotropy in the cosmic microwave background. The observations 
from cosmic background explorers differential radiometer have detected and measured 
cosmic microwave background anisotropics in different angular scales. These anisotropics 
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are supposed to hide in their fold the entire history of cosmic evolution dating back to the 
recombination era and are being considered as indicative of the geometry and the content 
of the universe. More about cosmic microwave background anisotropy is expected to be 
uncovered by the investigations of microwave anisotropy probe. There is widespread con- 
sensus among the cosmologists that cosmic microwave background anisotropics in small 
angular scales have the key to the formation of discrete structure. The theoretical ar- 
guments and recent experimental data that support the existence of an anisotropic 
phase that approaches an isotropic one leads to consider the models of universe with 
anisotropic background. A BI universe, being the straightforward generalization of the 
flat FRW universe, is one of the simplest models of an anisotropic universe that describes 
a homogeneous and spatially fiat universe. Unlike the FRW universe, which has the same 
scale factor for each of the three spatial directions, a BI universe has a different scale factor 
in each direction, thereby introducing an anisotropy to the system. It moreover has the 
agreeable property that near the singularity it behaves like a Kasner universe, even in the 
presence of matter, and consequently falls within the general analysis of the singularity 
given by Belinskii et al. Also in a universe filled with matter for = C < 1; 

it has been shown that any initial anisotropy in a BI universe quickly dies away and a 
BI universe eventually evolves into a FRW universe Since the present-day universe 

is surprisingly isotropic, this feature of the BI universe makes it a prime candidate for 
studying the possible effects of an anisotropy in the early universe on present-day ob- 
servations. In light of the importance mentioned above, several authors have studied BI 
universe from different aspects. In these papers we considered the nonlinear spinor field, 
as well as a system of interacting spinor and scalar fields. Beside the spinor fields, we also 
study the role of a A term in the evolution of the universe. 

As it is mentioned earlier, in our previous papers we limited our study within the scope 
of BI universe for it's simplicity and it's quick transition to a FRW universe in the process 
of evolution. But there are few other models, which describe an anisotropic space-time 
and generate particular interest among physicists [13] . The aim of this report is to study 
the self-consistent system of nonlinear spinor and anisotropic gravitational fields. As an 
anisotropic space-time we choose Bianchi type VI (BVI) model, since a suitable choice of 
it's parameters evokes few other Bianchi models including BI and FRW universe. It can 
be noted that unlike the BI universe, BVI space-time is an inhomogeneous one. Inclusion 
of inhomogeneity in the gravitational field significantly complicates the search for an exact 
solution to the system. In the Section 2, we write the equations for nonlinear spinor fields 
and the system of Einstein equations. In this section we also give their solutions in some 
general form, more precisely, we write the solutions in terms of a time-depended function 
that can be defined only given the concrete form of the spinor field nonlinearity. In 
Section 3, we give exact solutions to the equations for some special choices of spinor field 
nonlinearity and space-time inhomogeneity. Beside this, we also present some numerical 
solutions in graphical form. 



We shall investigate a self-consistent system of nonlinear spinor and Einstein gravita- 
tional fields. These two fields are to be codetermined by the following action: 



2. FUNDAMENTAL EQUATIONS AND GENERAL SOLUTIONS 




(2.1) 
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with 

L = Lgrav. "I~ ^spinor- (2-2) 

The gravitational part of the Lagrangian ()2.2|) is given by a Bianchi type- VI (BVl) space- 
time, while the spinor part is a usual spinor field Lagrangian with an arbitrary nonlinear 
term. 



A. Spinor field Lagrangian 

For a spinor field symmetry between \j/ and iff appears to demand that one should 
choose the symmetrized Lagrangian [14]. Keep it in mind we choose the spinor field 
Lagrangian with a nonlinear term in (j2.Hj) as follows: 



I 

^spinor ~ 



M\[f\{r + LN. (2.3) 



Here M is the spinor mass. The nonlinear term L/y describes the self- interaction of spinor 
field and can be presented as some arbitrary functions of invariants generated from the 
real bilinear forms of spinor field having the form: 

S = xffxi/ (scalar), (2.4a) 

P = i^ry^xif (pseudoscalar), (2.4b) 

= iwy^W) (vector), (2.4c) 

= (xfry^Y^xir) (pseudo vector), (2.4d) 

Q^"^ = (ij/c7^^i//) (antisymmetric tensor), (2.4e) 

where O^^ = (//2)[7^7^ — 7^7""]. Invariants, corresponding to the bilinear forms, are 

/ = 5^ (2.5a) 

J = /'^ (2.5b) 

Iv = Vf,v^ = {i[rr^yr)g^iv{vf¥). (2.5c) 

Ia = A^A^' = {xfr7^r^xir)g^,y{xfry^f¥). (2.5d) 

Iq = QiivQ^^ = (v/a^»g^a^v/3(V^cy"V)- (2.5e) 



According to the Pauli-Fierz theorem [1^ among the five invariants only / and J are 
independent as all other can be expressed by them: ly = —Ia = I + J and Iq = I — J. 
Therefore, we choose the nonlinear term L/^ = F(/,7), thus claiming that it describes the 
nonlinearity in the most general of its form. 



B. The gravitational field 

The gravitational part of the Lagrangian in ()2.2p has the form: 

_ R 

Lgia\. — (2.6) 
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Here R is the scalar curvature and K is the Einstein's gravitational constant. The gravi- 
tational field in our case is given by a Bianchi-type VI (BVI) metric: 

ds^ = dt^-a^e-^"''dx^-b^e^"'dy^-c^dz'^, (2.7) 

with a, Z?, c being the functions of time only. Here m, n are some arbitrary constants and 
the velocity of light is taken to be unity. It should be emphasized that the BVI metric 
models a Universe that is anisotropic and inhomogeneous. A suitable choice of m, n as well 
as the metric functions a, b, c in the BVI given by ()2.7|1 evokes the following Bianchi-type 
universes. Thus 

• for m = n the BVI metric transforms to a Bianchi-type V (BV) one, i.e., m = n, BVI 
=4* BV G open FRW with the line elements 

ds^ = dt^ - a^e-^"''dx^ - b^e^"''dy^ - c^dz^; (2.8) 



• for n = the BVI metric transforms to a Bianchi-type HI (Bill) one, i.e., n = 0, 
BVI =^ Bill with the line elements 

ds^ = dt^ - a^e-^'"' d^ - b^ dy^ - dz^; (2.9) 



• for m = n = the BVI metric transforms to a Bianchi-type I (BI) one, i.e., m = n = 0, 
BVI =^ BI with the line elements 

ds^ = dt^ - dx^ - b^ dy^ -c^dz^; (2.10) 



for m = n = and equal scale factor in all three directions the BVI metric transforms 
to a Friedmann- Robertson- Walker (FRW) universe, i.e., m = n = and a = b = c, 
BVI =^ FRW with the line elements. 

ds^ = dt^-a^{dx^ + dy^ + dz^). (2.11) 



Let us write the nontrivial components of Ricci and Riemann tensors as well as Christof- 
fel symbols of the BVI metric. 

The nontrivial Christoffel symbols of the BVI metric read 

Toi = a/a, = b/b, T^g = c/c, 

yl - _^ r2 - « r3 - "^p-^'^z y3 _ 2mz 

1 31 — m, I32 — n, 1 11 — -^e , 1 22 — 

The nontrivial components of Riemann tensor are 

pOl _ ^ p02 _ ^ p03 _ ^ 
^01 - ^02--^, ^03---, 

12 _ mn db ,3 _ cd .3 _ be 
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The nontrivial components of the Ricci tensor are 

n f d b , ,c 

= [m n [m — n]- 

\ a b c 



1 /d db dc rrP' — mn 

R\ = - - + -T + ^ 

\a ab ac 
9 /b db be rP' — mn 

R? = - T + -T + 



^ \b ab be e^ 

;c dc be m^ + n^ 

m = - - + — + 



c ' a c b c 

To investigate the existence of singularity (singular point), one has to study the invari- 
ant characteristics of space-time. As we know, in general relativity, these invariants are 
composed from the curvature tensor and the metric one. Though, in 4D Riemann space 
there are 14 independent invariants, it is sufficient to study only three of them, namely 
the scalar curvature l\ =R, h =R^vR^'^, and the Kretschmann scalar = Rap^vR*^^^^ ■ 
From the Riemann and Ricci tensors written above one can easily write find 

2 

Ii=R = — \f — dbc — abc — dbc — ab{m^ — rnn + ri^) /c], T = abc, (2-12) 

T 

whereas I2 °« 1/(tc)^ and oc 1/(tc)^. Thus we see that at any space-time point, where 
T = the invariants Ii,l2,h become infinity, hence the space-time becomes singular at 
this point. 



C. Field equations 

The field equations for the spinor and gravitational fields can be obtained from vari- 
ational principle. Variation of the Lagrangian ()2.3p with respect to field functions Yiv) 
gives the nonlinear spinor field equations: 

iY^'V^\j/-M\j/+^\j/ + M)^\j/ = (2.13a) 
/V^V/7^+Mv/-^V/-/^V/7^ = (2.13b) 

where ^ = 25"^/ and ^ = 2PFj. 

Varying ()2.3p with respect to metric function {gfj,v) we find the Einstein's field equation 

Rli-^dl^R = KTy^, (2.14) 

where Ry is the Ricci tensor, R is the Ricci scalar, and Ty is the energy-momentum 
tensor of the spinor field. In our case, where space-time is given by a BVI metric ()2.7p . 
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b c be 




= 


-+-+--- 




b c be 




c d c d 


2 

m 




- + -H 




= 


c a c a 






d b db 


mn 




- + + + 






a b ab 







the equations for the metric functions a, b, c read 

(2.15a) 
(2.15b) 
(2.15c) 

db be cd m^ — mn + n^ _o 

-T + T- + 2 = '^^0' 2.15d 

a b b c c a 

m--n--(m-n)- = kT^. (2.15e) 
a b c 

Here over dots denote differentiation with respect to time (t). The energy-momentum 
tensor of the material field has the form: 

7? = y'''' (^n^v¥+¥rv^^i¥-^^i¥rv¥-^vn¥^ -s^i^p. (2.16) 

Here L^p is the spinor field Lagrangian, which on account of spinor field Eqns. fl2.13p 
takes the form: 

Lsp = -^S-'^P + F. (2.17) 

In the expressions above denotes the covariant derivative of spinor, having the form 
0: 

y^V=dfiyr-r^X{f (2.18) 

where is spin connection. The spin affine connection matrices F,, (x ) are uniquely 
determined up to an additive multiple of the unit matrix by the Eqn. [l7l | 

dfiYv-K^ra-T^rv + rv^fi = o, (2.19) 

with the solution 

TmW = \8po{x){d,eieP-rP^,)ff. (2.20) 
For the metric element ()2.7|1 it gives 

To = 0, 

Fi = -\dff-m-ff]e-'"' 
2 '- c 

T2 - ^[bff + n^ff]e"\ 



It is easy to show that 



where we define 



x = abc. (2.21) 
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The Dirac matrices Y^{x) of curved space-time are connected with those of Minkowski 
one as follows: 

yO = f, Y^ = fe"''/a, f = f/be"\ Y^ = f/c 

with 



f 



I 
-/ 



f 



where Oi are the Pauli matrices: 

.1 _ / 1 
1 



a' 
-a' 



-/ 
/ 



, y^ = f 



-/ 
-/ 



1 

-1 



Note that the y and the o matrices obey the following properties: 

ffJ + fJf = 2ri'J, /,j = 0,l,2,3 
ff + ff = 0^ {ff=I, / = 0,1,2,3 
o^o^ = 5jk + iejMo\ J, fc, / = 1 , 2, 3 

where 77,^ = {1, — 1, — 1, — 1} is the diagonal matrix, 5jk is the Kronekar symbol and £jif_i 
is the totally antisymmetric matrix with £123 = +1. Let us consider the spinors to be 
functions of t and z only, such that 

V/(^,^)=v(Oe*^ W{t.z)=v{t)e-''' 
Inserting ()2.22|) into ()2.1H|1 for the nonlinear spinor field we find 



(2.22) 



m — n k\ ^ ^ r 



V + — V 

2t 



f 



Ic 
m — n 



+ i-]v'f' -i^v-'^vf = 0. 
2c c/ 



(2.23a) 
(2.23b) 



Here we define ^ — M — Let us introduce a new function 

Uj{t) = V^Vj{t). 

Then for the components of the NLSF from ()2.23p one obtains 



Ml + i^ui — 
U2 + i^U2 + 
U3 — Z$M3 — 
U4 — /4>M4 + 



m — n k 

— /-+^ 

2c c 

m — n k ^ - 

2c c 
m — n .k 

2c c 

m — n k 

— /-+^ 

2c c 



"3 

Ul 
U2 



0, 
0, 
0, 
0. 



(2.24a) 
(2.24b) 
(2.24c) 
(2.24d) 



Using the spinor field Eqns. ()2.13|1 and ()2.23|1 it can be shown that the bilinear spinor 
forms, defined by ()2.4|) . i.e., 

S= ij/i// = vv, p = ixfry^xir = ivfv, = xfrfy^x^/ = vfy^v, 
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obey the following system of equations: 



So-2^Ql^-2m'^Q = 0, (2.25a) 
/b-2^<2o^-24>A° = 0, (2.25b) 



^A^ + 24)Po + 2^5o = 0, (2.25c) 
aI —AI = 0, (2.25d) 

Vq-^Vq = 0, (2.25e) 
Vi - + 2*^30 _ 2^g2i ^ 0^ (2.25f ) 

Qf + 2-^SQ-2^Vi = 0, (2.25g) 

Qf + 2-^PQ + 2mi = 0, (2.25h) 

where we denote Fq = xF . Combining these equations together and taking the first integral 
one gets 

{So? + {Po? + {Alf - {Alf - {V^f + (y3)2 + (230)2 ^ ^g2i )2 ^ ^ ^ Const (2.26) 

Before dealing with the Einstein Eqns. ()2.15|) let us go back to ()2.24|) . From the first 
and the third equations of the system ()2.24j) one finds 

"13 = (^ - Q)u\^ - 2i^un + (^ + 2), (2.27) 

where, we denote mo = ui/u-^ and Q = [m — n — 2ik]/2c. The Eqn. (j2.27|) is a Riccati 
one ^3 with variable coefficients. A transformation (3| 

vi3=exp(- j{<^-Q)undt), (2.28) 

leads the general Riccati Eqn. ()2.27j) to the second order linear one, namely, 

(^ - 2) V13 + [2/4>(^ _ 2) _ + 2] + (^^ - e)2(?f + 2)vi3 = 0. (2.29) 

Some times it is easier to solve a linear second order differential equation than a first 
order nonlinear one. But we also give a general solution to ()2.27|) . For this purpose we 
rewrite ()2.27j) in the form 

VV13 = (^-2)w?3e-2'^*W* + (^ + 2)e2'7*W*^ (2.30) 

where we set M13 = Wi3exp[— 2// ^{t)dt]. It is an inhomogeneous nonlinear differential 
equation. The solution for the homogeneous part of ()2.3m) . i.e., 

VV13 = (^ - e)wl3e-2'7*(0^' (2.31) 
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reads 



Wl3 = - 



(2.32) 



with C being some arbitrary constant. Then the general solution to the inhomogeneous 
Eqn. ()2.3()j) can be presented as 



Wi3 = - 



1 -1 



(2.33) 



with time dependent parameter C{t) to be determined from 



C 



j{^-Q)e-^'S^^')d'dt + C{t) 



' + Q)e 



2iJ(i>{t)dt 



[2M) 



Thus given the concrete nonlinear term in the Lagrangian and solution of the Ein- 
stein equations, one finds the relation between ui and M3 {u2 and U4 as well), hence the 
components of the spinor field. 

Now we study the Einstein Eqns. fl2.15|) . In doing so, we write the components of the 
energy-momentum tensor, which in our case read 



mS-F + -V\ 
c 



= -kV'. 

Let us demand the energy-momentum tensor to be conserved, i.e.. 



Taking into account that is a function of t only, from ()2.36p we find 



k . -I km — fici 
4)5o - 5f Po + -Vi 



c c 
m — n. 



-Vk 



0, 
0. 



(2.35a) 
(2.35b) 
(2.35c) 
(2.35d) 

(2.36) 

(2.37a) 
(2.37b) 



As one can easily verify, the Eqns. ()2.37|) are consistent with those of ()2.25p . 
Let us go back to the Eqns. ()2.15j) . In view of ()2.35j) . from ()2.15e|l one obtains 



(^-^ ^exp[-Kk j V^dt], ^ = const. 
Subtracting ()2.15a|) from ()2.15b|) we find 



(2.38) 



d 
dt 



2 2 
m — n 



(2.39) 
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Analogically, subtraction of ()2.15a|) from ()2.15c|) and ()2.15b|) from ()2.15c|) give 



and 



d 
dt 

d 

dt 



d ( /a 
dt\ \c 



mn + m 



Kk 



Kk T 

— y^T, 

c 



(2.40) 



:2.41) 



respectively. It can be shown that, in view of (IT^ and (ITip the Eqns. (ESHI), 
and ()2.41|) are interchangeable. Thus in ()2.38|) we have inter relation between the metric 
functions, whereas two more equations for defining any two of the three. Taking into 
account that T = abc, from (12.381) we can write a and b in terms of c, such that 



x"c"'-^''^exp[-Kk / V^dt] 



1 / (m+n) 



and 



/{^exp[-Kk / V^dt]) 



l/{in+n) 



In view of (EH, and (jT^ from one finds 



Kk 



(m + n) 



[2A2) 



(2.43) 



(2.44) 



In getting ()2.44p we employ only four out of five Einstein equations, leaving ()2.15d|) 
unused. 

On the other hand, adding (I2.15aj) . (I2.15bj) , (I2.15c|) and (I2.15djl . multiplied by 3 we 
get the equation for T, which in view of ()2.35|1 takes the form 



K 



+ - [3{MS + ^5 + - 2F) + 2-y3] . 



:2.45) 



Thus we are left with two equations, namely ()2.44|) and ()2.45|) for two unknowns c and T. 
These two equations can be combined to get 



^ + — = — y3 + ^ — + M5 + ^5 + ^/'-2F . 

c X c c 2 



:2.46) 



Thus we come to the Eqn. (j2.46p where all the equations in hand, both spinor and 
gravitational, are employed. Moreover, both c and T are inter-related. Assuming c as a 
function of T (vice versa) and giving the concrete form of spinor field nonlinearity one 
finds the solution of ()2.4(jj) exactly what we do in the section to follow. 



3. ANALYSIS OF THE RESULT 

In the previous section we derived the fundamental equations for nonlinear spinor fields 
and metric functions. Comparing the equation with those of in a BI universe (see e.g., jof) 
we conclude that introduction of inhomogeneity both in gravitational (through m and n) 
and spinor (through k) fields significantly complicates the whole picture. In what follows. 
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we try to write the solutions in a more explicit form for some special choice of spinor field 
nonlinearity and field and space-time inhomogeneity. 

Let us first consider the Eqn. ()2.46p . As one sees, there are two unknown functions in 
this equation, namely, c and T. Though, they are inter-related, there is no equation giving 
this relation explicitly. So, as a first step, we have to make some assumption relating c 
and T. In what follows, we consider few models giving explicit relation between c and T 
and study them for different types of nonlinear spinor terms. 

Case I. Let us assume that 

c = T. (3.1) 

Under this assumption in view of T = abc, we should have a= l/b. Indeed, from ()2.42|) 
and we find 

r f n 1 l/iin+n) 

a = x"'-"^cxp[-Kk / V^dt] 



and 



exp[-Kk I V'^dt]) 



l/(m+;j) 



With regard to ()3.1|) . from ()2.46|) we obtain 



f Kk rn^ + n^ K r ^ ^ -, 
- = —V^ + — ^ + -\MS + ^S + WP-2F]. 



(3.2) 
(3.3) 

(3.4) 



Let us now study ()3.4j) for some special choice of spinor field nonlinearity. 
Linear spinor field 

To begin with we consider the linear case setting F{I,J) = 0. It immediately leads to 
^ = and ^ = 0. Eqn. ()3.4|) now takes the form 



X Kk rn^ + n^ K 

- = —V^ + ^ + -MS. 

XX x^ 2 

73 



(3.5) 



As one sees, to solve ()3.5p . we have to find V and S first. From ()2.25|) for the linear spinor 
field we have 



So-2-Ql^ 
Po-2'^qI'-2MA^q 



m — n 



■AI + 2MPq 



■ m — n r, 



m — n 



V^ + 2MQ'^ 



30 



Q^^ + 2-SQ-2MVi 



0, 
0, 
0, 
0, 
0, 
0, 
0, 
0, 



(3.6a) 
(3.6b) 
(3.6c) 
(3.6d) 
(3.6e) 
(3.6f) 
(3.6g) 
(3.6h) 
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with the first integrals 

iSo? + {V^? + {Ql'f-{V^? = 0, (3.7a) 
(^o)2 + (A«)2 + (2^1)2 _(a3)2 = 0. (3.7b) 

Thus we see that even in case of linear spinor field with k we cannot write or 5 
explicitly. In order to express S or P, hence the massive term or spinor field nonlinearity, 
in terms of T, we now consider the spinor field to be to be space independent setting k = 0. 
From ()2.25j) in this case one obtains, 



Co 



(3. 



with Co being the integration constant. 

The Eqn. fl3.4|) in this case takes the form 



T 2 

The solution of ()3.9p can be written in quadrature as 

dz 



+ ^MCo. (3.9) 



^2(m2 + n2) Inr + kMCq t + E 



const. (3.10) 



For the solution to be meaningful, the integrand in ()3.10p should be positive. It means 
that for the T to have an initial value close to zero, one has to set small value for m and 
n, while the constant E should be large enough. 

Components of spinor field can be obtained from ()2.32|) . In the case considered, = 0, 
Q= {m-n)/lx and 4> = M. 

Nonlinear spinor field with k = Q. 

Let us now consider the nonlinear spatially independent spinor field. We first choose 
the nonlinear term being a function of / = ^2 only, followed by a massless spinor field with 
the nonlinear term to be a function of J = P^. 

If the nonlinear spinor term is given as F = F{I) = XS^ , where A is the (self) coupling 
constant, then in view of 5 = Cq/t for T we find 

^=— + jMC„+ ZJ" - ("1) 

with the solution in quadrature 

dx 



2(m2 + n2)ln T + kMCq t - kAcJ t2-'7 +E 



t. (3.12) 



As one sees, inclusion of nonlinear term sets additional restriction on the smallness of the 
initial value of T. Components of spinor field, as in linear case, can be obtained from 
()2.32|) . In the case considered, ?f = and Q — {m — n) /Ic. 

For the massless spinor field, if the nonlinear tern is chosen as F = F{J) = XP^ , from 
(IT^ for P we find 

P = Dq/x. (3.13) 
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FIG. 1: Perspective view of T for power of nonlinearity T] = 4 and the integration constant 
E = I. Here "I" stands for nonlinear case with F = F{I), "J" for massless spinor field with the 
nonhnearity type F = F{J) and "L" for linear spinor field. 



Equation for T then takes the form 

+ \'J (3.14) 

with the solution in quadrature 



m^ + n^ fcA(7]-2)Dj 



2(m2 + n2)ln T - fcAD^ T^-^ + E 



t. (3.15) 



Components of spinor field can be obtained from ()2.32j) . In the case considered, Q — 
(m — n)/2c and 4> = 0. 

Here we illustrate some numerical results obtained for cases considered above. The 
parameters of the equations are taken as follows: for spatial inhomogeneity parameters 
we set m = 0.0002 and n = 0.0001, whereas Einstein's gravitational constant K is taken to 
be unity. For the nonlinear spinor field we choose A = 0.5. 

As one sees from Fig. ^ the presence of the massive term accelerates the expansion 
process. For the given parameters in case of linear spinor field, the initial value of T can 
be chosen very small (here we set Tq — l.e — 6), whereas for the nonlinear cases there exists 
lower boundary, which for the set of parameters are Tq = 0.57 and Tq = 1.0, for "I" and 
"J" cases, respectively. Note that for different value of r\ the corresponding pictures are 
similar to Fig. ^ with the permissible initial value of T increasing with that of r\ . Note 
that the figures to follow are plotted for the nonlinear term being a function of invariant 
/, i.e., F=F(/). 

As is seen from Fig. |2l for a positive r/, the greater is the value of E the faster the 
corresponding T expands, though the permissible initial value of T for a greater E is 
less than that for a smaller one. In this particular case, we have Tq = 0.77, 0.57, 0.46 
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FIG. 2: Perspective view of T for T] =4 and different value of integration constant E. 



for E = 0.1, 1.0, 2.0, respectively. Here we also note that given the negative value of 7] 
it is possible to obtain oscillatory mode of expansion. We do exactly that for the case 
considered below. 

Case II. Let us now consider the case setting 



c = \fx. 

This leads to following expressions for a and h 



a 



J V^dt] 



1 / (m+«) 



and 



r'V2 



exp[-Kk / V^dt]) 



1 / (m+n) 



Under this assumption from ()2.46p we get 

f = 2KkV^y/T + 2{m^ + n^) + K [MS + ^S + ^P- 2F] x. 



(3.16) 

(3.17) 
(3.18) 

(3.19) 



This equation can be solved exactly as in previous case if we set = and choose the spinor 
field nonlinearity as F = F{I) or in case of massless spinor field F = F{J) or F = F{I± J). 

For the reason that will be given afterwards, we consider the nonlinear spinor field in 
a Bianchi type V (BV) universe setting m = n in the corresponding equations. To begin 
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with we write the equations for bilinear spinor forms. Setting m = n in ()2.25|) one finds 



T/3 





= 0, 


/'o-24>A|] 


= 0, 




= 0, 


^0 


= 0, 




= 0, 




= 0, 




= 0, 




= 0, 



with the following relations between spinor bilinear forms 



0\2 



{VoY + iQ'or + iQo 



30\2 



52, 
53, 
^4, 



(3.20a) 
(3.20b) 
(3.20c) 
(3.20d) 
(3.20e) 
(3.20f) 
(3.20g) 
(3.20h) 

(3.21a) 
(3.21b) 
(3.21c) 
(3.21d) 



where 5, are the constant of integration. 

Let us now go back to Einstein's equations. Eqn. (|2.15ep in this case takes the form 



d b 
a b 



m 



(3.22) 



Unlike the BVI, where the corresponding equation, i.e., ()2.15e|) inter-connects all the three 
metric functions a, b, c, Eqn. ()3.22|) relates only a and b between them: 

a = ^exp[-{Kk/m) Jv^dt]b. (3.23) 

Recalling T = abc in view of ()3.23j) we can now express a and b in terms of T and c: 



b = ^-^/^^/T/dexp[{Kk/2m) jv^dt]. 
In view of (jT^ . (jT^ and the fact that = 0, from 



d_ 
dt 



dt 



{<)}] 



2m^ 



one obtains 



- = 3-- 



, m 



On the other hand ()2.45|1 in this case has the form 

,2 



T 2 c 



(3.24) 
(3.25) 

(3.26) 
(3.27) 

(3.28) 



16 



Bijan Saha 



Combining ()3.27|) and ()3.28p we obtain 

^ + — = —V^ + 2^ + -\MS + ^S + ':^P-2F]. 3.29 

c T c 3c c^ 2 

Thus we see that a strait-forward insertion of m = n into ()2.46j) does not lead to ()3.29j) . 
since Eqn. ()2.15e|) for different Bianchi type space-time gives different relations between 
the metric functions. Here we simply note that for a Bill metric, where n = 0, Eqn. 
()2.15e|l relates a with c, whereas for a BI universe, as well as for FRW universe there is 
no such equation. Note that, though in a BV space-time where m = n, many equations 
in question become significantly simpler, it is not enough to write the solutions explicitly, 
smce V\ S and P are still undefined explicitly. As in previous case, we again consider 
only time dependent spinor filed setting k = 0. It will give us enough ground to solve both 
spinor and gravitational field equations explicitly. 

Before studying the Eqn. ()3.29|) in details, we go back to nonlinear spinor field equa- 
tions. With m = n and k = for the spinor field we immediately find 

iii -f-z4>Mi -^M3 = 0, (3.30a) 

U2 + i^U2-'^U4 = 0, (3.30b) 

M3-i4>M3+^Ml = 0, (3.30c) 

U4-i^U4 + '^U2 = 0. (3.30d) 

As in BVI space-time we consider the nonlinear term to be F = F{I), or for massless 
spinor field F = F{J) or F = F{I± J). The spinor field Eqn. ()3.30|1 completely coincides 
with those for a BI metric. So in what follows we simply write the corresponding results 
without any details. A detailed analysis of these results can be found in 

Thus, for the nonlinear terrn_in the Lagrangian given a.s F = F{I), the components of 
the spinor field take the form 





= (Ci/v^)exp[- 


-ifiM 


-9)dt], 


(3.31a) 


V2{t) 


= (C2/v^)exp[- 


-ijiM 


-^)dt], 


(3.31b) 


V3{t) 


= (C3/v^)exp[? 


j{M- 


^)dt], 


(3.31c) 


¥4it) 


= (C4/v/T)exp[« 


j{M- 


9)dt]. 


(3.31d) 



Here Ci, C2, C3, C4 are the integration constants, such that 

-f Cf — C3 — C4 = Co, 



with Co = St. 

In case, the nonlinear term is given by F = F(7), the components of the spinor field 
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lift — 

V'l — 




¥2 = 


-^{D2e''' + iD4e 
V ^ 


¥3 = 


-^{iDie''' +D3e 


¥4 = 


-^{iD2e''' +D4e 



have the form 

■^), (3.32a) 

(3.32b) 
(3.32c) 

^). (3.32d) 

Here o — J ^dt, and the integration constants D,- obey 

2{dI + dI-DI-DI)=Do, 

with Dq to be determined from P = Dq/x. 

Thus we see that in the cases considered here the spinor bilinear forms are inverse 
proportional to T, i.e., S = Cq/t and P = Dq/t. 

Let us now go back to (|3.29p . As one sees, for k = 0, the assumption t = c makes no 
sense, since in this case the metric functions a and b turns out to be constant. So as was 
mentioned earlier, we consider the case with c = \fx. Under this assumption from ()3.29p 
we get 

f = 4m2 + fc[M5 + ^5 + ^/'-2F]T. (3.33) 

We consider the case with F = XS^ . Taking into account that 5 = Cq/t and ^ = 0, from 
()3.33|) one derives 

x = 4m^ + K [MCq + XC^ {ri - 2) t^"'^] , (3.34) 
with the solution in quadrature 

t. (3.35) 



(8m2 + 2kMCo) t - 2A KCjj t2-'7 + E 

Note that, for a linear spinor field and for the massless spinor field with F = XP^ one has 
to put A = and M = 0, respectively, into ()3.35|1 . It should be noted that for a positive 
constant in case of linear spinor field T may have even a trivial initial value. 

Finally we give numerical solutions in graphical form for case II. It should me mentioned 
that for a positive r\ we have the picture that we obtained in case I with the expansion 
being faster in this case. So we concentrate our study for the negative r\ only. In this 
case we leave the value of inhomogeneity parameters and Einstein constant unaltered. 

As one sees from Fig. 01 a negative 77 gives rise to the oscillatory mode of evolution. For 
the given set of parameters, for E < 0.001 we have the solution that is almost periodical, 
i.e., in this case the model begins to expand as soon as it reaches the singular point 
corresponding T = 0, whereas for greater E these phases are separated by a longer time 
interval. 

In figure Fig. |3]we plot the evolution of T for different rj with E being fixed. It should 
be mentioned that unlike the case with positive 77 in this case solutions are available only 
for even rj. Finally we like to state that choosing a negative 77 we come to the analogical 
conclusion for case I. 
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FIG. 3: Perspective view of T for a negative 7] with different E. In this particular case we set 
7] = —4. 




o 2 -4. e s -io 12 1-4 ie 

Time 

FIG. 4: Perspective view of T for £ = 0.001 with different value of T). 

4. CONCLUSION 

Self-consistent system of nonlinear spinor field and gravitational one, modeled by a 
Bianchi type VI (BVI) space-time, is studied. Exact solutions to spinor and gravitational 
field equations are obtained for some special choice of spinor field nonlinearity and space- 
time inhomogeneity. It is shown that if the nonlinear spinor term is chosen to be a 
function of the invariants I = or J = P^, with a negative power, the model provides 
oscillatory mode of expansion. It is shown that though a suitable choice of m and n in 
a BVI metric evokes other Bianchi models, namely, BV, Bill and BI, the solutions of 
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Einstein equations in these universes cannot be obtained by simply setting m and n in the 
corresponding solutions obtained in a BVI universe, since in different models the metric 
functions are inter-related differently. Indeed, it follows from the Eqn. ()2.15e|) . 




that for a BVI model the metric functions a, b, c are connected with each other by ()2.38|) . 
whereas, for BV universe ()4.H1 gives relation between a and b by ()H.23j) and for a Bill 
space-time it connects a and c. For BI or FRW models Eqn. ()4.1|) does not exist. 



Yu.P. Rybakov, B. Saha and G.N. Shikin, PFU Reports: Physics 2, (2), 61 (1994). 
Yu.P. Rybakov, B. Saha and G.N. Shikin, Commun. Theor. Phys. 3, 199 (1994). 
B. Saha and G.N. Shikin, J. Math. Phys. 38, 5305 (1997). 

R. Alvarado, Yu.P. Rybakov, B. Saha and G.N. Shikin, JINR Preprint E2-95-16, 11 p. 
(1995), Commun. Theor. Phys. 4, (2), 247 (1995), gr-qc/9603035 

R. Alvarado, Yu.P. Rybakov, B. Saha and G.N. Shikin, Izvestia VUZov: Fiz. 38,(7) 53 
(1995) [Russ. Phys. J. 38, 700 (1995)]. 

B. Saha, and G.N. Shikin, Gen. Relativ. Gravit. 29, 1099 (1997). 
Bijan Saha, and G.N. Shikin, PFU Reports: Physics 8, (1), 17 (2000); |g^/ 0102059[ 
Bijan Saha, Mod. Phys. Lett. A 16, 1287 (2001); gr-qc/0009002j 
Bijan Saha, Phys. Rev. D 64, 123501 (2001); gr-qc/0107013, 

C. W. Misner, Astrophys. J. 151, 431 (1968). 

V.A. Belinskii, I.M. Khalatnikov and E.M. Lifshitz, Adv. Phys. 19, 525 (1970). 
K.C. Jacobs, Astrophys. J. 153, 661 (1968). 

T. Singh and A.K. Agrawal, Int. J. Theor. Phys. 32, 1041 (1993). 
T.W.B. Kibble, J. Math. Phys. 2, 212 (1961) 

V.B. Berestetski, E.M. Lifshitz and L.P. Pitaevski, Quantum Electrodynamics (Nauka, 
Moscow, 1989). 

V.A. Zhelnorovich, Spinor theory and its apphcation in physics and mechanics (Nauka, 
Moscow, 1982). 

D. Brill and J. Wheeler, Rev. Mod. Phys. 29, 465 (1957). 

E. Kamke, Differentialgleichungen losungsmethoden und losungen (Akademische Verlags- 
gesellschaft, Leipzig, 1957). 

[19] V.F. Zaitsev and A.D. Polyanin, A Handbook on Nonhnear Differential Equations (Nauka, 
Moscow, 1993). 



